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Abstract 

When impurity and phonon scattering coexist, the Boltzmann equation has been solved accu- 
rately for nonlinear electron transport in a quantum wire. Based on the calculated non-equilibrium 
distribution of electrons in momentum space, the scattering effects on both the non-differential (for 
a fixed dc field) and differential (for a fixed temperature) mobilities of electrons as functions of 
temperature and dc field were demonstrated. The non-differential mobility of electrons is switched 
from a linearly increasing function of temperature to a parabolic-like temperature dependence as 
the quantum wire is tuned from an impurity-dominated system to a phonon-dominated one [see 
T. Fang, et ai, Phys. Rev. B 78, 205403 (2008)]. In addition, a maximum has been obtained in 
the dc-field dependence of the differential mobility of electrons. The low-field differential mobility 
is dominated by the impurity scattering, whereas the high-field differential mobility is limited by 
the phonon scattering [see M. Hauser, et al, Semicond. Sci. Technol. 9, 951 (1994)]. Once a 
quantum wire is dominated by quasi-elastic scattering, the peak of the momentum-space distribu- 
tion function becomes sharpened and both tails of the equilibrium electron distribution centered 
at the Fermi edges are raised by the dc field after a redistribution of the electrons is fulfilled in 
a symmetric way in the low-field regime. If a quantum wire is dominated by inelastic scattering, 
on the other hand, the peak of the momentum-space distribution function is unchanged while 
both shoulders centered at the Fermi edges shift leftward correspondingly with increasing dc field 
through an asymmetric redistribution of the electrons even in low- field regime [see C. Wirner, et 
al, Phys. Rev. Lett. 70, 2609 (1993)]. 
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I. INTRODUCTION 



A clear understanding of the electron transport processes in both semiconductor quan- 
tum wells and quantum wires is of importance for the development of high- mobility quantum 
electronic devices. Several studies have been connected with the low-temperature electron 
conduction in quantum wires within the linear-transport regime. ^"^^ The resonant effects of 
both magnetic field and electron-electron interaction on thermoelectric power in quantum 
wires were investigated. I0'i3,i4 presence of electron-electron interaction varying along 

quantum wires with a broken translational symmetry, the linear temperature dependence of 
resistivity was predicted. When a voltage is applied between spht gates, the backscatter- 
ing of electrons at one of two rough edges inside a conduction channel was found enhanced 
experimentally, although the phenomenon has not been fully understood so far in a perpen- 
dicular magnetic field due to combined effects of the Lorentz force and a static potential 
across split gates. For the effect of electron-phonon interaction on conduction carriers in 
quantum wires, non-Markovian Langevin-like equations were derived for calculating carrier 
temperature and a current-voltage instability was predicted. ^^'^^ Moreover, the Boltzmann 
transport equation was solved for linear transport of electrons in single and tunnel-coupled 
double quantum wires with a full inclusion of electron-phonon scattering on a microscopic 
level. ^1 

In recent times, the linear magneto-quantum transport properties of electrons in single 
and tunnel-coupled double quantum wires have been extensively investigated by Lyo and 
Huang ^'^^'^^ based on an accurate solution of the linearized Boltzmann equation in which 
the effects of both elastic and inelastic scattering as well as the effect of electron-electron 
scattering on the conductance of electrons had been treated microscopically beyond the well- 
known relaxation-time approximation Furthermore, the nonlinear electron transport 
under a strong dc field in a single quantum wire was also explored by using the coupled force 
balance and microscopic scattering equations in which only the elastic scattering effect 
on the electron mobility had been computed. The replacement of the energy-balance equa- 
tion for the relative scattering motion of electrons by the microscopic scattering equation 
allows one to incorporate the asymmetric electron distribution in momentum space. How- 
ever, the microscopic scattering equation still employs an adiabatic-type approximation ^^"^^ 
to the coupling between the center-of-mass and relative motions of electrons. The aim of 
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present paper is to develop an efficient numerical method for solving accurately the Boltz- 
mann transport equation which includes inelastic scattering beyond the linear-transport 
regime and further applying it to quantum -wire systems as a first step. 

Prom the current study, we find that the coupled force balance and microscopic scattering 
equations are still not accurate enough to describe the electron distribution as a function of 
momentum in the high-field regime. Although the calculated currents as functions of either 
dc-field or temperature can be compared directly with the experimentally measured ones, 
the momentum dependence of the electron distribution in the nonlinear-transport regime 
cannot be extracted reliably from the measured current on the microscopic level. Therefore, 
the accurate calculation of the steady-state non-equilibrium electron distribution based on 
the Boltzmann equation becomes necessary for this purpose and is expected to provide 
a fair evaluation of the accuracy of the coupled force-balance and microscopic scattering 
equations in the high-field limit. In this paper, we modify the previously derived formalism 
by Huang, Lyo and Gumbs^^ for the Bloch oscillations of electrons within a quantum-dot 
superlattice to investigate the nonlinear transport of electrons in a single-subband quantum 
wire. Here, for completeness, we will include some of the key equations from our previous 
studies The destruction of phase coherence by electron-phonon scattering in disordered 
conductors has been investigated. However, we shall confine our attention to high-quahty 
quantum wires where the localization length is expected to be longer than the sample length. 
Compared with a strong electron-phonon interaction, the electron-electron interaction here 
is negligible for a quantum-wire system with low electron densities. As long as the subband 
spacing is much larger than the average thermal energy as well as the Fermi energy of con- 
duction electrons, the single-subband model employed here is justifiable. ^ This condition is 
easily achieved for a quantum- wire system with a strong lateral confinement and low electron 
densities at moderate to low lattice temperatures. As is well known, electron-electron inter- 
action yields the Tomonaga-Liittinger liquid effect in clean quantum-wire systems. Indeed, 
a conductance increase by electron-phonon interaction in quantum wires was predicted for 
this case. However, this effect will not be considered here for diffusive electron transport. 

Our main results in this paper include the comparisons of quasi-elastic with inelastic 
scattering effects on both the temperature T and the dc-field jFg dependence of electron mo- 
bility, as well as on the wave number /c-dependence of the electron extreme non-equilibrium 
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distribution function. For the electron mobility as a function of T, its linear increase with T 
in an impurity-dominated system "^^ switches to a parabolic- like T dependence in a phonon- 
dominated system with a maximum occurring at a characteristic temperature at which the 
strengths of quasi-elastic and inelastic scattering of electrons become comparable with each 
other. In the presence of scattering, we have found that the electron mobility increases 
(decreases) with in the low (high) field regime, resulting in a maximum in the mobility 
as a function of !Fq. In addition, when the quantum wire radius R is decreased, the char- 
acteristic temperature is also reduced due to the enhanced phonon scattering strength, and 
the maximum in the J^o-dependent mobility shifts down to a lower value. Furthermore, 
a well-known 1/T reduction of the mobility in the high-T limit ■^^ is obtained when the elec- 
tron density is lowered due to more accessible phonons contributing to a scattering process, 
and the mobility is increased (decreased) for low (high) values because of the reduced 
impurity scattering (the enhanced nonlinear phonon scattering) of electrons, respectively. 

For electron distribution as a function oi k/kp (k-p is the Fermi wave number), the peak at 
/c = is unchanged while both shoulders centered at k = ±kp shift leftward with increasing 
To in a quantum-wire system dominated by inelastic scattering. For a system dominated by 
quasi-elastic scattering, on the other hand, the peak of the distribution function at A; = 
becomes sharpened. At the same time, both tails of the equilibrium electron distribution 
in the region > kp are lifted up by the quasi-elastic scattering of electrons. Our calcu- 
lated asymmetric (symmetric) distribution functions of electrons for low-field regime in the 
quantum-wire systems dominated by inelastic (quasi-elastic) scattering cannot be described 
by a quasi-equilibrium Fermi function employed in the energy balance equation with an 
effective electron temperature 

The paper is organized as follows. In Sec. II, we introduce our model for single-subband 
quantum wire systems, derive the general formalism for nonlinear electron transport, and 
solve the Boltzmann equation with the full inclusion of scattering by both impurities and 
acoustic phonons. In Sec. Ill, numerical results are presented for comparison between the 
effects for quasi-elastic and inelastic scattering of electrons on the temperature and dc- 
field dependence of mobility, as well as for a comparison of the quasi-elastic and inelastic 
scattering effects on the electron distributions in momentum space. The conclusions drawn 
from these calculated results are briefly summarized in Sec. IV. 
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II. MODEL AND THEORY 



In this section, we describe a microscopic approach for the Boltzmann equation in a quasi- 
one-dimensional quantum wire. Scattering by both impurities and acoustic phonons is fully 
taken into account. We assume a strong parabolic potential for the transverse confinement 
of the electrons. 

We begin by considering a degenerate electron gas moving in a quantum wire in the y 
direction, as shown in Fig. 1. For this system with a strong transverse confining potential, 
only the lowest subband will be occupied by electrons in the low-density and low-temperature 
regime, and therefore, a single subband model will be adequate. The original form of the 
Boltzmann equation for a time-dependent distribution function f{k, t) of electrons in a 
quantum wire is given by 



+ Y. (^+(*. *') (Ka + 1) /(*'. *) [1 - /(*. ')] - *) [1 - /(*'. 

+ V-ik, k') f{k', t) [1 - f{k, t)] - + 1) f{k, t) [1 - f{k\ t)]}) , , (1) 

where jFg denotes the applied dc electric field, k is the electron wave vector along the y 
direction, q is the three-dimensional wave vector of phonons, \ = ^ {\ = t) labels the 
longitudinal (transverse) acoustic phonon modes, and nq\ — NQ{u!qx) is the Bose distribution 
function for the equilibrium phonons with frequency cUgx in the mode of q and A. In Eq. (1), 
we have defined the electron-phonon inelastic scattering rate as ^^'^^''^'^ 

V±{k, k') = — |Vfc,fe'p5(£fc - Ek' ± hujgx) k',k±qy , (2) 
and the electron- impurity elastic scattering rate as ^'^^ 

27r 

Ve{k,k') = —\Uk,k'\'S{ek-ek'), (3) 

where Ek — ti^k"^ /Im* is the electron kinetic energy with m* the electron effective mass, Vk,k' 
is the interaction between electrons and phonons, Uk,k' is the interaction between electrons 
and impurities, and the so-called [/-scattering process has been neglected due to low 
densities and low temperatures of the system. 



We define f(k, t) — fj^^ +gk(t) with a dynamical non-equilibrium part gk(t) for the total 
distribution function, where fjf^ = fo{£k) is the Fermi function for equilibrium electrons. 
We note that the terms of zeroth order in gk in the scattering parts of Eq. (1) cancel out due 
to detailed balance. After incorporating these considerations, we can algebraically simplify 
Eq. (1) into 



dgk{t)_ 
dt 



dsk 



+ ^ + [9-k{t) - 9k{t)] k') 

(0) 



- 9k{t) E [P+ik, k') (n,, + /I?) + V-{K k') (n,A + 1-/1' 

fc',q,A 

+ gk^{t)[V^{k,k')-V-{k,k')] 

+ E 9k'{t) k') (n,x + 1 - /f ) + V-{K k') (n,A + /f ) 

fe',q,A 



(4) 



where Vk = {l/h)dek/dk = hk/m* is the group velocity of electrons with the wave number 
k and the /c'-sum for k' e —k is carried out near the resonance at k' = —k. 

The transport relaxation rate from electron-impurity scattering equals ^ 



- = 2 E K{k,k') = l\Uk,-k\'Viek) . 



(5) 



k'e-k 



Here, the factor 2 arises from the 2kF {kp is the Fermi wave vector at zero temperature) 
scattering across the Fermi sea and 'D{e) is the total density-of-states given by 



V{e) 



1 2m* 



(6) 



nh \ Ek ' 

where Lq is the length of the quantum wire. For electron-phonon scattering, we define 
scattering rates: 



>v±(^, k') = E r±{k, k') . 

q,A 



(7) 



Making use of Eq. (2) and writing hujgx = ax^J {k — k'Y + q\ from the Debye model, we 
have^^ 



(8) 
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where 9{x) is the unit step function, S is the sample cross-sectional area, ax — hcx, cx is 
the sound velocity for the phonon in the A mode, \qy\ — \k' — k\, 



Q± = V^' + = ;r V (^'^ - ^k'T - o^xik - k'Y , (9) 



? = (10) 

(y.\ 

and the form factor in Eq. (8) takes the following form 

= exp(-?ii?72) = exp {-[{sk - Sk'f - o?x{k - k'f]R^/2al} (11) 

for a transverse parabolic confining potential within the a;2;-plane with the radius R for 
the cylindrical confinement model. The quantity ^{qi_) arises from the consequence of 
momentum conservation for g^. A smaller value of R in Eq. (11) implies a stronger electron- 
phonon scattering in Eq. (8). 

Assuming axial symmetry, the strength of the electron-phonon interaction |VA(q)P in 
Eq. (8) is given by ^° 



m%)fpciv ^ ■ I „ , , (12) 

where V = SLq is the sample volume, p is the mass density, D is the deformation-potential 
coefficient, h\i is the piezoelectric constant, e{qy) is the Thomas-Fermi dielectric function, 
and is the structure factor for acoustic phonons due to anisotropic electron-phonon 
coupling. The quantum- wire structure factors of phonons in Eqs. (12) are calculated to be 

^'^^ - 2(gi + g^)3 ' %i + g,^)^ ■ ^ ^ 

We further note from Eq. (8) that only one of W±(/i;, A;') is nonzero, i.e., either 
>V+(/c, A;') > for e^/ > £fe or yV-{k^ k') > for Sk' < Sk- This directly yields the re- 
sult 



yV±{k, k') = 0(±£fe, T Sk) J2 >Va(A:, k') , , (14) 

A 
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where W\{k, k') — Wa(^', k) is given by 



m{k, k') = ^e \{eu - sw? - al{k - k'f \ |l^,(q)p ^ . (15) 

11 ^ ^ ax 

Using the above scattering rates, we rewrite Eq. (4) in the concise form 



- 9k{t) E k') [nu^y + + ^0 + 1 - ff 

k' 

+ gu'it) [W+{Kk')-W-{Kk')]} 

+ E ^^'(0 K(A^, k') (rife,,, + 1 - /f) + A;') (n^,,, + /f) 



(16) 



where n^^fe/ = A^od^fc ~ independent of A = By introducing the notations 

fk — fk'^ fk — ^ ~ fk'^ obtain the discrete form for the total inelastic scattering 
rate 



= ^ E >^±(^' + f^) ' 



i',± 



with A; = jSk and A;' = j'5k {5k — > 0), and 



We can also write 



>V|(t) = ^ 5^ E W [>^+(^' k') - W.{k, k')] , 



and 



(17) 



(18) 



(19) 



dt 



eJ^o dgk{t) 



dsk h dk 
9k{t) [W,- + Wlit)] + ^ [^_,(i) - gk{t)] + ^SkY: 9k'{t) W 



3,3 



(20) 



In Eq. (20), the quantity fiffeWj is the only nonlinear term in g^. 

Equation (20) is equivalent to the following matrix equation for 1 < j < 
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dt 



N 



i'=i 



(21) 



where k = [j - {N + l)/2] 5k, k' = \j' - {N + l)/2] 5k, 5k = 2k^^/{N - 1) with AT > 1 an 
odd integer, and < /Cmax- The elements for the vector b in Eq. (21) are 



bj — eJ^oVj ^ 



(22) 



and the elements for the matrix a(t) in Eq. (21) are 



h5k 



1-5 



'j, (JV+l)/2 



2Ti 



-5. 



<j+j',N+l 



1-5 



'j,{N+l)/2 



2Ti 



Cjj/ — — 5k Wi. ,v . 



27r 



(23) 



By assuming that the dc electric field is turned on at time t — 0, the differen- 
tial equation (21) can be solved in combination with the initial condition 51^(0) = for 
j — 1, 2, • • • , N. In addition, using the so-called 3-point central-difference formula, which 
transforms the differential term (ej^o/h) dgk{t)/dk in Eq. (20) into a finite-difference term, 
we can write Cjji in Eq. (23) in a simple form of 



(24) 



Unfortunately, the quantities {gi, g^, • ■ ■ , gN} are not linearly independent of each other, 
which implies that the matrix a.{t) in Eq. (21) is a singular one. In fact, the condition for 
the particle-number conservation requires that 



N 



12 9At) = o. 

Combined with the fact that Sk — £-fe and Vk — —V-k, this yields for 1 < j < 



(25) 



N-l 

9N{t) = - E 9At) , (26) 

where j — N is the right- most k-space point. As a result of Eq. (26), we can renormalize 
the singular {N x A'")-matrix a.{t) in Eq. (21) into a regular [{N -1) x {N - l)]-matrix a'{t) 
through the following relation: 
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a'jj,{t)^aj,f{t)-aj,N{t) , 
where j,j' = 1, 2, • • • , iV — 1. Consequently, Eq. (21) is renormalized into 



(27) 



-^-b'^-j:/.At)9'At), (28) 

where the vectors g'{t) and b' are the same as the vectors g(t) and b, respectively, without 
the last element with j = N. In a similar way, we correspondingly rewrite Eq. (19), by 
excluding the element with j — N, as 



>V|(t) = fSkY: g'A^) {W+(j, f) - W_(j, f) - [W+(j, N) - >V_(j, N)]} . (29) 
j'=l 



Once the non-equilibrium part, gj{t), of the total electron distribution function has been 
solved from Eq. (28), the transient drift velocity, Vc{t), of the system can be calculated 
according to 



Vc{t) 



N 



(0) 



-1 



(30) 



The steady-state drift velocity of electrons is given by Vc{t) in the limit as i — > oo. 
Furthermore, the steady-state conduction current is given by Ic — eniDfd, where riio is the 
linear density of electrons in a quantum wire. The standard non-differential mobility for a 
fixed dc field is given by the ratio /Xnd = v^/ J^q, while the so-called differential mobility of 
electrons for a fixed temperature is defined through — dv^/dTg. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

For convenience of notation, we define tj — at the energy Sj — Sk in Eq. (5): 




where 70 is the value of l/r^ at Sj — E-p, and Ep is the chemical potential /Iq of electrons 
at T = OK. For a fixed electron temperature T and a linear density niD, A*o is determined 
from the following relation: 
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""'^ TT expiisj - iio)/kBT] + 1 ■ ^^^^ 

The Fermi wave vector is equal to kp — 7miD/2 and the Fermi energy is 
at T = K. We chose GaAs as the host material in carrying out our numerical calculations 
and made use of the following parameters: the sound velocities q = 5.14 x 10^ cm/sec, 
Ct = 3.04 X 10^ cm/sec, the mass density p = 5.3g/cm^, the piezoelectric constant hu = 
1.2 X lO^V/cm, and the deformation-potential coefficient D = — 9.3eV. For the elastic 
scattering, we chose 70 = 1 x 10^^ s~^. The other parameters, such as R, Tq, T and riiD, 
will be directly indicated in the figure captions. 



A. Mobility 

For linear electron transport, the non-differential mobility /i^d only depends on the sam- 
ple geometrical size, impurity and phonon scattering strengths, electron density niD, and 
temperature T but not on the dc field jFg. On the other hand, the differential mobility /j.^ de- 
pends on J^Q for nonlinear electron transport. In this case, the electron mobility approaches 
one constant at the low- field limit, while saturating at another constant in the high-field 
regime. Therefore, the J^o-dependent electron mobility implies a nonlinear transport pro- 
cess of electrons in a system. 

Figure 2 compares the quasi-elastic (very large values for R, blue dashed curves) with 
the inelastic (R — 30 A, solid red curves) scattering effects for riiD = 2 x 10^ cm~^ on the T 
dependence of non-differential mobility [in (a) when — 40 V/cm] defined by //nd — ^d/^o, 
as well as on the dependence of differential mobihty [in (b) with T = 20K] defined by 
Hd = dvd/dJ-'o. For the T dependence in Fig. 2(a), fi^d exhibits a simple linear relation when 
only the quasi-elastic scattering (i.e., with a very large value for R) exists. However, after 
the inelastic electron scattering has been introduced to the system, /ind decreases strongly, 
and the straight line is replaced by a parabola-like curve with a maximum occurring at a 
characteristic temperature. The characteristic temperature corresponds to the case when 
the strengths of both quasi-elastic and inelastic scattering become comparable to each other. 
Similar thermal behavior was predicted by Fang, et ai, [30] for the linear electron transport 
in graphene nanoribbons. On the other hand, for /id as a function of J-'q in Fig. 2(b), we 
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find dfXd/dTo > for the small values of To but dji^/dJ^Q < for high when either the 
quasi-elastic or the inelastic scattering of electrons is present. A strong reduction in at 
low is found for phonon scattering, however, this difference becomes much smaller at 
high To. 

In Fig. 3, we display the inelastic scattering effects at = 2 x lO^cm"^ on the T 
dependence of //nd [in (a) with — 40 V/cm] and Tq dependence of ji^ [in (b) with T = 20K] 
for electrons in quantum wires with i? = 10 A (solid red curves) and i? = 30 A (blue 
dashed curves). When R is decreased from 30 A to 10 A in Fig. 3(a), i.e. the phonon- 
scattering strength is increased, the system dominated by the quasi-elastic-scattering with 
d/ind/dT > for small values of T switches to one dominated by the inelastic-scattering with 
d/jind/dT < because of the decrease of the characteristic temperature with decreasing R. 
For the To dependence of //d in Fig. 3(b), the dc- field value corresponding to the maximum 
in //d is shifted downwards as the radius R of the quantum wire is decreased to 10 A, in 
combination with a strong reduction in /i^ at all the field values shown in this figure. 

From Figs. 2(b) and 3(b) we find that //d, as a function of To, initially increases but 
eventually decreases with To when the inelastic scattering effect is not very strong. In 
this case, the frictional force acted on electrons can be reduced ^^'^^'^^ as long as the lattice 
temperature T, as well as the impurity density, are low and the electron drift velocity 
Vd is less than the sound velocity Cx. The reduced frictional force to moving electrons is 
associated with the softened acoustic-phonon mode due to the Doppler shift. However, 
/Id of hot electrons decreases significantly once the nonlinear phonon scattering becomes 
dominant with increasing To- 

Figure 4 presents the results when inelastic scattering effects are present for i? = 10 
A on the T dependence of /x^d [in (a) with To — 40 V/cm] and the J'-'o-dependence of /id 
[in (b) with T = 20K] for the quantum wires with niD — 2 x lO^cm"^ (solid red curves) 
and riiD = 1 x lO^cm"^ (dashed blue curves). When nio decreases from 2 x lO^cm"^ 
to 1 X lO^cm"^ in Fig. 4(a), more acoustic phonons in the system become accessible for 
the inelastic scattering of electrons since ujqx ~ 2cxkF oc niD in this case. As a result of 
this, the well-known 1/T reduction of fi^d at the high-temperature limit ^'^ is seen. At the 
same time, /Xnd is enhanced for low T but reduced for high T with decreasing rin). The 
enhancement of /i^d at low T is due to the reduced I/tj oc ^/EF oc riiD in Eq. (31) for the 
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quasi-elastic scattering of electrons, while the reduction of //nd at high T directly comes from 
increased phonon scattering proportional to T once ^ hcuqx. A small enhancement in 
/Xd for low J^o and a large reduction in /i^ for high J^q are also found in Fig. 4(b) as tt-id 
is reduced to 1 x 10^ cm^^. In this case, however, the de-field value corresponding to the 
maximum in is not shifted. The little enhancement of /^d at low JFq is attributed to the 
decreased quasi-elastic scattering of electrons, while the significant reduction of fid at high 
J^o is a consequence of the increased nonlinear phonon scattering [W|(t) in Eq. (29)] which 
directly contributes to the J^o dependence of /i^ for nonlinear electron transport. In this case, 
the low-field mobility is dominated by impurity scattering, while the high-field mobility is 
limited by phonon scattering. The increase in the high-field mobility with electron density 
was experimentally observed previously by Hauser, et ai, [33] in a quantum wire. 

B. Distribution Function 

Prom Eq. (30) we are aware that the steady-state drift velocity is uniquely determined 
by the non-equilibrium part g'^it) of the total electron distribution function f{k, t) in the 
limit t ^ oo. Prom a physical point of view, this implies that there must be a microscopic 
origin for all the observed behaviors for nonlinear electron transport in Figs. 2(b), 3(b) and 
4(b), and this origin is believed to be related to the field- induced redistribution of conduction 
electrons in /c-space. 

In Fig. 5, we show the inelastic scattering effects for R — 10 A, T = 10 K and tt-id = 
2 X 10^ cm~^ in both the non-equilibrium part of electron distribution function gk [in (a)] 
and the total electron distribution function fk [in (b)] as functions of a scaled electron wave 
number k/kp for To — 10 V/cm (black solid curves), 20 V/cm (red dashed curves), 40 V/cm 
(green dash-dotted curves), and 60 V/cm (blue dash-dot-dotted curves). For the sake of 
comparison, we also plot the Fermi function f^^^ (thick black dashed curve) in Fig. 5(b). In 
Fig. 5(a), gife > {gk < 0) represents the population (depopulation) of fc-state electrons in the 
quantum wire. When JFq is gradually increased from 10 V/cm up to 60 V/cm, the inelastic 
phonon scattering in the system moves more and more electrons from the right Fermi edge 
at k/kp — +1 to a high-energy states peaked at k/kp — —2 after passing through the left 
Fermi edge at k/kp — —1. However, the electrons at the lowest k — energy state are 
left largely untouched. This significant adjustment of conduction electron distribution in 
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k-space can be also seen from the total distribution function fk in Fig. 5(b), where the peak 
at A; = is unchanged while both shoulders centered at k/kp — ±1 are shifted to the left 
with increasing J^q. This leads to an asymmetric distribution function with respect to A; = 
even in the low-field regime, a non-quasi-equilibrium electron distribution. A similar Fermi- 
edge shifting in a hot-electron distribution with an increasing dc field was experimentally 
observed by Wirnel, et al, [34] in GaAs/AlGaAs heterostructures. 

Figure 6 exhibits the effects of quasi-elastic electron scattering at T = 10 K and 
riiD = 2 X lO^cm"^ on gt [in (a)] and fk [in (b)] as functions of k/kp for JFq = lOV/cm 
(black solid curves), 20V/cm (red dashed curves), 40V/cm (green dash-dotted curves), 
and 60V/cm (blue dash-dot-dotted curves). For clarity, we also include fj^^^ (thick black 
dashed curve) in Fig. 6(b). As J^q increases from 10 V/cm up to 60 V/cm in Fig. 6(a), all the 
conduction electrons in the low-energy k states between the left and right Fermi edges are 
depopulated and moved to high-energy k states outside the region with \k/kp\ < 1. However, 
the depopulation of electrons at = is found to be smaller than those at the two Fermi 
edges. This specific redistribution of electrons in A';-space by quasi-elastic scattering is also 
reflected in Fig. 6(b) for fk, where the peak at A; = is sharpened and the electrons around 
the peak are moved outward to the high-energy tail states of the equilibrium distribution 
function fjf^ in A;-space. However, moving electrons outward by quasi-elastic scattering in 
both directions in A;-space creates a somewhat symmetric distribution of electrons in the 
low-field regime with a sharpened peak and two raised tails at the same time, in strong 
contrast to the result for inelastic scattering in Fig. 5(b). Although more and more high- 
energy states become populated by moving electrons from low-energy states through both 
quasi-elastic scattering and inelastic scattering of electrons, moved electrons in high-energy 
states only have negative group velocities, i.e. dsk/dk < 0, for the inelastic scattering case. 



IV. CONCLUDING REMARKS 



The accurate solution of the Boltzmann equation for a quantum wire has been obtained 
with the full inclusion of both impurity and phonon scattering for conduction electrons on the 
microscopic level. Using this solution, both the non-differential and differential mobilities 
of the electrons as functions of temperature and dc field, as well as the non-equilibrium 
distribution of electrons in momentum space, were analyzed. 
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As a result of our studies, we have found that the non-differential mobility of electrons 
in a quantum wire switches from a linear increase with T in an impurity-dominated system 
to a parabolic-like T dependence in a phonon-dominated system. The characteristic tem- 
perature, at which the inelastic and quasi-elastic scattering strengths become comparable 
to each other, corresponds to a maximum in the non-differential mobility as a function of 
temperature, and it shifts with the quantum-wire radius. We also find that the differential 
mobility of electrons has a maximum as a function of dc field, and the field corresponding 
to this maximum shifts as the quantum- wire radius is varied. The low-field mobility is dom- 
inated by impurity scattering, while the high-field mobility is limited by phonon scattering. 

For a quantum-wire system dominated by quasi-elastic scattering, the peak in the dis- 
tribution function at A; = becomes sharpened, and both tails of the equilibrium electron 
distribution centered at the Fermi edges are lifted up by the dc field. On the other hand, for 
a quantum wire dominated by inelastic scattering, the peak at /c = is unchanged, while 
both shoulders centered at the Fermi edges shift to the left correspondingly with increasing 
dc field. Therefore, the quasi-elastic scattering in a quantum wire redistributes the conduc- 
tion electrons in a symmetric way in momentum space in the low-field regime. However, 
the inelastic scattering of electrons redistributes them in an asymmetric way even in the 
low-field regime. 

Acknowledgment s 

The authors would like to thank the Air Force Office of Scientific Research (AFOSR) for 
its support. 



[1] J. Lee and H. N. Spector, J. Appl. Phys. 54, 3921 (1983). 
[2] G. Fishman, Phys. Rev. B 34, 2394 (1986). 

[3] H. Sasaki, T. Noda, K. Hirakawa, M. Tanaka, and T. Matsusue, Appl. Phys. Lett. 51, 1934 
(1987). 

[4] S. Das Sarma and X. C. Xie, Phys. Rev. B 35, 9875 (1987). 



15 



[5] K. F. Berggren, G. Roos, and H. van Houten, Phys. Rev. B 37, 10118 (1988). 

[6] H. Akira and T. Ando, Phys. Rev. B 43, 11676 (1991). 

[7] B. Tanatar and A. Gold, Phys. Rev. B 52, 1996 (1995). 

[8] S. K. Lyo and D. Huang, Phys. Rev. B 64, 115320 (2001). 

[9] T. Sugaya, J. P. Bird, M. Ogura, Y. Sugiyama, D. K. Ferry, and K. Y. Yang, Appl. Phys. 
Lett. 80, 434 (2002). 

[10] S. K. Lyo and D. H. Huang, Phys. Rev. B 66, 155307 (2002). 

[11] S. K. Lyo and D. H. Huang, Phys. Rev. B 68, 115317 (2003). 

[12] S. K. Lyo and D. H. Huang, J. Phys.: Condensed Matter 16, 3379 (2004). 

[13] A. M. Lunde, K. Flensberg and L. Glazman, Phys. Rev. Lett. 97, 256802 (2006). 

[14] A. M. Lunde, K. Flensberg and L. Glazman, Phys. Rev. B 75, 245418 (2007). 

[15] J. Rech and K. A. Marveev, Phys. Rev. Lett. 100, 066407 (2008). 

[16] D. Hartmann, L. Worschech and A. Forchel, Phys. Rev. B 78, 113306 (2008). 

[17] L. G. Mourokh, J. Phys.: Condens. Matter 9, 8489 (1997). 

[18] L. G. Mourokh, Phys. Rev. B 57, 6297 (1998). 

[19] S. K. Lyo and D. H. Huang, Phys. Rev. B 73, 205336 (2006). 

[20] J. M. Ziman, Principles of the Theory of Solids, pp. 212 -215, (Cambridge Press, Second 
edition, 1972). 

[21] Z.-G. Wang, D. Suqing, and X.-G. Zhao, Phys. Rev. B 69, 035305 (2004). 

[22] X.-G. Zhao, G. A. Georgakis, and Q. Niu, Phys. Rev. B 56, 3976 (1997). 

[23] C. S. Ting, S. C. Ying, and J. J. Quinn, Phys. Rev. B 14, 4439 (1976). 

[24] D. H. Huang and G. Gumbs, Phys. Rev. B 80, 033411 (2009). 

[25] X. L. Lei and C. S. Ting, Phys. Rev. B 32, 1112 (1985). 

[26] X. L. Lei, Phys. Rev. B 77, 033203 (2008). 

[27] D. H. Huang, S. K. Lyo and G. Gumbs, Phys. Rev. B 79, 155308 (2009). 

[28] J. Rammer and A. Schmid, Phys. Rev. B 34, 1352 (1986). 

[29] T. Brandes and A. Kawabata, Phys. Rev. B 54, 4444 (1996). 

[30] T. Fang, A. Konar, H. Xing, and D. Jena, Phys. Rev. B 78, 205403 (2008). 

[31] D. H. Huang, T. Apostolova, P. M. Alsing, and D. A. Cardimona, Phys. Rev. B 69, 075214 
(2004). 

[32] D. H. Huang, P. M. Alsing, T. Apostolova, and D. A. Cardimona, Phys. Rev. B 71, 195205 

16 



(2005). 

[33] M. Hauser, E. Gornik, C. Wirner, M. Baur, G. Bohm, and G. Weimann, Semicond. Sci. 
Technol. 9, 951 (1994). 

[34] C. Wirner, C. Kiencr, W. Boxleitner, M. Witzany, E. Gornik, P. Vogl, G. Bohm, and G. 
Weimann, Phys. Rev. Lett. 70, 2609 (1993). 



17 



CD 
O 

o 



Quantum Wire 





Curren^ 






>• 

y 




X 




FIG. 1: (Color online) Schematic illustration for a quasi-one-dimensional quantum wire driven by 
an electric field along the y (wire) direction, where the confinements in the z and x directions are 
given by a cylindrical parabolic potential. The current flowing through the channel between source 
and drain electrodes is indicated by a horizontal (red) arrow. 
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FIG. 2: (Color online) Calculated non-differential mobility ^nd = Vd/^o in (a) ^ function of 
temperature T and differential mobility /id = dv^/dJ-Q in (b) as a function of dc-field J-q for 
inelastic scattering with i? = 30 A, (solid red curves) and quasi-elastic scattering (blue dashed 
curves) at niD = 2 x 10^ cm^^. In (a), we take J-q = 40 V/cm, while we choose T = 20K in (b). 
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FIG. 3: (Color online) Numerical values of /ind in (a) as a function of T and fi^ in (b) as a function 
of JFq for two values of i? = 10 A, (solid red curves) and i? = 30 A (blue dashed curves) at 
niD = 2 X 10^ cm^^. In (a), we take J^o = 40 V/cm, while we choose T = 20K in (b). 
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FIG. 4: (Color online) Calculated /ind in (a) as a function of T and /id in (b) as a function of J-q 
for different electron densities nio = 2 x 10^ cm^^ (solid red curves) and nio = 1 x 10^ cm^^ (blue 
dashed curves) with i? = 10 A. In (a), we take J^o = 40 V/cm, while we choose T = 20K in (b). 
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FIG. 5: (Color online) Steady-state gk in (a) and fk in (b) as functions ol k/k-p for various dc-fields: 
J^o = lOV/cm (black solid curves), 20V/cm (red dashed curves), 40V/cm (green dash-dotted 
curves), and 60V/cm (blue dash-dot-dotted curves). Here, we take R = 10 A, T = 10 K and 
niD = 2 X 10^ cm~^. For the sake of comparison, we also show fj^^ (thick black dashed curve) in 
(b). 
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FIG. 6: (Color online) Plots of steady-state gk in (a) and fk in (b) as functions of k/kp for various 
dc-fields: J^o = lOV/cm (black solid curves), 20V/cm (red dashed curves), 40V/cm (green dash- 
dotted curves), and 60V/cm (blue dash-dot-dotted curves). Here, we take T = 10 K, tiid = 
2 X 10^ cm^^ and include the quasi-elastic scattering. For a clear comparison, we also display fj^^ 
(thick black dashed curve) in (b). 
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